A circular dielectric grating for vertical extraction of single quantum dot emission 



O 

(N 

Oh! 



C/5 1 

O ' 

H— > . 

Ph. 

O ; 

o ■ 

• i-H , 

43 ■ 

Oh' 



> 

(N 
O 

o 



X 



M. Davanco, 1 ' 2 ^ M. T. Rakher, 1 D. Schuh, 3 A. Badolato, 4 and K. Srinivasan 1 

1 Center for Nanoscale Science and Technology, National Institute of Standards and Technology, Gaithersburg, MD 20899, USA 

2 Maryland NanoCenter, University of Maryland, College Park, MD 
3 Institute for Experimental and Applied Physics, 
University of Regensburg, D-93053 Regensburg, Germany 
^Department of Physics and Astronomy, University of Rochester, Rochester, New York 14-627, USA 

(Dated: January 18, 2013) 

We demonstrate a nanostructure composed of partially etched annular trenches in a suspended GaAs 
membrane, designed for efficient and moderately broadband (~ 5 nm) emission extraction from 
single InAs quantum dots. Simulations indicate that a dipole embedded in the nanostructure center 
radiates upwards into free space with a nearly Gaussian far field, allowing a collection efficiency 
> 80 % with a high numerical aperture (NA=0.7) optic, and with m 12 x Purcell radiative rate 
enhancement. Fabricated devices exhibit a ~ 10 % photon collection efficiency with a NA=0.42 
objective, a 20 x improvement over quantum dots in unpatterned GaAs. A fourfold exciton lifetime 
reduction indicates moderate Purcell enhancement. 



Efficient extraction of single photons emitted by indi- 
vidual semiconductor epitaxial quantum dots (QDs) is a 
necessity for many applications in spectroscopy and clas- 
sical and quantum information processing [l| . As epitax- 
ially grown QDs are embedded in semiconductor mate- 
rial, total internal reflection of the emitted light at the 
semiconductor-air interface and radiation divergence can 
typically lead to < 1 % collection efficiencies even with 
high numerical aperture (NA) optics. Photonic struc- 
tures such as micropillar cavities can provide both QD ra- 
diative rate Purcell enhancement and a far-field radiation 
pattern that can be effectively collected 2], but require 
precise spectral tuning of the cavity resonance to the 
QD emission line. In contrast, vertically-oriented etched 
nanowires 0) are spectrally broadband structures that 
have recently been shown to provide large free space col- 
lection efficiencies, albeit without Purcell enhancement 
and with an involved fabrication process. Broadband op- 
eration not only relaxes the spectral alignment require- 
ment, which may impose strict constraints in fabrication 
tolerances, but is also a necessity in spectroscopic appli- 
cations in which simultaneous detection of various spec- 
trally separate transitions is desired. Here, we present an 
approach for efficient free space extraction of QD emis- 
sion using a suspended circular grating. This structure 
requires a simple nanofabrication procedure, and sup- 
ports a relatively broad (few nanometer) optical reso- 
nance with a directional, nearly Gaussian far-field, which 
allows efficient free space photon collection. Simulations 
predict a collection efficiency of w 53 % (80 %) into a 
NA=0.42 (0.7) optic. In fabricated devices, we report a 
« 10 % single QD photoluminescence (PL) collection effi- 
ciency into a NA=0.42 objective, a«20x improvement 
compared to QDs in unpatterned bulk GaAs. A four- 
fold reduction in QD lifetime is also observed, indicating 
moderate radiative rate enhancement. 
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FIG. 1: (a) Top, (b) angled, and (c) cross-sectional SEM 
images of suspended circular dielectric grating structure. 



Our nanostructure (Fig. [1} consists of a circular di- 
electric grating with radial period A that surrounds a 
central circular region of radius 2A, produced on a sus- 
pended GaAs slab of thickness t — 190 nm. The GaAs 
slab supports single TE and TM polarized modes (elec- 
tric or magnetic field parallel to the slab, respectively). 
The grating is composed of ten partially etched circu- 
lar trenches of width w and depth d, with t/2 < d < t. 
Quantum dots are grown at half the GaAs slab thickness 
(z = 0), and located randomly in the xy plane. This 
'bullseye' geometry favors extraction of emission from 
QDs in the central circular region. It is based on (linear) 
high-contrast second-order Bragg gratings recently intro- 
duced [4] for light extraction from planar waveguides. 
While similar circular geometries have been employed for 
enhanced light extraction from light emittingdiodes 
and for demonstrating annular Bragg lasers [6|, here we 
show an application in QD single photon extraction. 

The design process consisted of a series of finite dif- 
ference time domain simulations that maximized vertical 
light extraction near the expected QD s-shell emission 
(^QD ~ 940 nm), by varying A, t, and w. The struc- 
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FIG. 2: Electric field intensity in the (a) xy and (b) xz planes 
(log scale), (c) Far- field polar plot for the cavity mode with 
A = 360 nm. (d) Collected power (P co i) as a function of 
varying NA, normalized by the upwards (Pz+) and total (Ptot) 
emitted powers, (e) Calculated vertically extracted power 
as a function of wavelength, normalized to the homogeneous 
medium electric dipole power PHom for d = 0.70t. Continuous 
lines: upwards (+z) extraction; dotted: downwards (—2). (f) 
Experimental PL spectra for high QD density devices. 



tures were excited with a horizontally oriented electric 
dipole at the bullseye center (x = 0,y = 0), representing 
an optimally placed QD. Total radiated power, steady- 
state upwards emission, and electromagnetic fields were 
then recorded at several wavelengths. The grating period 
A was initially chosen to satisfy the second-order Bragg 
condition, A = Xqd/it-te, to allow for efficient vertical 
light extraction (ute is the GaAs slab TE mode effective 
index). The dipole orientation was assumed to be aligned 
along the xy plane, exciting only TE slab waves. Start- 
ing values for trench width and depth were w — 100 nm 
and d — 0.5/;, deemed to be easily fabricated. Vertical 
light scattering at the gratings is partial, so that second- 
order Bragg reflections towards the center lead to verti- 
cally leaky cavity resonances as shown in Figs. [5Ja) and 
(b). The large index contrast at the trenches leads to 
strong reflections and out-of-slab-plane scattering at the 
semiconductor-air interfaces, evident in the strong field 
concentration at the bullseye center in Fig. [3Ja) and the 
fast field decay within the first couple of trenches from 
the center (Fig. HJb)). Large differences in propagation 
constants in the semiconductor and air produce signif- 
icant resonance spectral shifts with small variations in 



trench width. Trench depth (d) has a strong influence 
on the quality factor (Q) and vertical light extraction, as 
incomplete spatial overlap between a trench and an inci- 
dent slab-bound wave leads to both coupling to radiating 
waves and lower modal reflectivity. Preferential upwards 
vertical extraction results from the grating asymmetry, 
and is optimized through the trench depth Q- We note 
that in addition to the mode shown in Fig. [TJ the cav- 
ity supports resonances which can be excited by dipoles 
offset from the bullseye center. Coupling to these reso- 
nances can lead to modified spontaneous emission rates 
and collection efficiencies Q. 

Figure[2je) shows simulated, upwards (continuous) and 
downwards (dotted) vertically extracted power as a func- 
tion of wavelength for structures with A = 350 nm, 
360 nm, and 370 nm, w = 110 nm, and d 0.70i. 
All curves are normalized to the homogeneous medium 
electric dipole power, PHom- Trench parameters reflect 
a trade-off in cavity Q and vertical light extraction, as 
discussed above. It is apparent that for each A, an 
« 5 nm wide resonance exists, with preferential upwards 
(+z) light extraction. The upwards extracted power is 
~ 10 x PHom, an indication of Purcell radiation rate en- 
hancement due to the cavity 0. Indeed, for the A = 360 
nm structure, on which we now focus, the enhancement 
F p at the maximum extraction wavelength (A c = 948.9 
nm) is F p = Ptot/PHom = 11-0, where P to t is the to- 
tal radiated power in all directions. This resonance has 
Q = 200, and its effective mode volume, calculated from 
the field distribution, is V e g = 1.29(A c /n) 3 (n is the GaAs 
refractive index) [7[ . The value for F p predicted by Q and 
V e f[ is wll.8, and is consistent with the value determined 
above by the dipole radiation simulations. Note that, 
given the modal field distribution in Fig. (Ha) , the mod- 
ified dipole emission rate depends strongly on its spatial 
location, being maximal at the bullseye center. 

The steady-state fields at a surface just above the 
GaAs slab were used to calculate the far-field pattern in 
Fig. (He), which shows that the emission is nearly Gaus- 
sian with a small divergence angle. To better quantify 
this, we calculate the power P co i collected by an optic of 
varying NA. Figure [U[d) shows the fractions of the up- 
wards emitted (P z +) and total (Ptot) powers collected as 
a function of the collection optic acceptance angle. For 
NA=0.42 (24.8° acceptance angle), » 60 % of the up- 
wards emitted power (or » 53 % of the total emission) 
can be collected. For NA> 0.7, or an acceptance angle 
> 44.4°, collection superior to 80 % of the total emission 
can be achieved. We note that our suspended grating 
approach limits radiation into the substrate without the 
need to oxidize the AlGaAs, bond the grating to a low 
index layer or utilize a deeply etched geometry @, Q. 

Gratings were fabricated in a t — 190 nm GaAs layer 
containing a single layer of InAs QDs (density gradi- 
ent from > 100 /im~ 2 to /jm~ 2 along the (Off) di- 
rection) on top of a I /im thick Alo.6Gao.4As sacrifi- 
cial layer Fabrication steps included electron-beam 
lithography, plasma dry etching, and wet chemical etch- 
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FIG. 3: (a) PL spectrum from a low QD density A = 360 nm 
device, for various pump powers, (b) Temperature evolution 
of spectrum in (a) (25 nW pump), (c) Temperature evolution 
of excitonic energies. Continuous lines are fits. 



tra. Saturated photon rates (collected with a NA=0.42 
objective) from Xi and X2 were at least 20 times higher 
than from typical QDs embedded in unpatterned GaAs, 
as shown in Fig.QJa). Assuming 100 % QD quantum ef- 
ficiency, we estimate a collection efficiency of « 10 % 
is achieved with the bullseye pattern Q. A lifetime 
measurement of Xi after a « 300 pm bandpass filter 
(Fig.rjjb)) exhibits a multi-exponential decay with a fast 
lifetime of 360 ps, limited by the ~ 600 ps timing jitter 
of the detectors. For comparison, the lifetime of a single 
QD inside of a suspended GaAs waveguide [l3[ (dotted 
in Fig. HJb)), for which no radiative rate modification is 
expected, was « 1.5 ns. This suggests F p > 4. Note 
that since the pump in Fig. @|a) is pulsed with a 20 ns 
repetition period, significantly longer than the lifetime, 
the increase in detected counts relative to unpatterned 
GaAs is solely due to enhanced photon extraction and 
collection into the objective. 



ing. The plasma dry etch was optimized so that the GaAs 
would be partially etched to a desired depth in the grat- 
ing region (Fig. [IJb)), and fully etched over the curved 
rectangles just outside the grating region (Fig. QJa)), 
which were used in the wet etching step to undercut and 
suspend the device. 

Testing was done in a liquid He flow cryostat at « 8 K. 
Figure^f) shows PL spectra of three devices with a high 
QD density and A =350hm, 360 nm, and 370 nm and 
d/t > 0.7, for pulsed pumping at a 780 nm wavelength 
(above the GaAs bandgap). The spectra closely resem- 
ble the theoretical curves of Fig.^e), with three, 5 nm 
wide peaks spaced by « 20 nm. Deviations are likely due 
to differences in geometry and refractive index between 
simulated and fabricated structures. These results vali- 
dated our simulations, and served to calibrate the fabri- 
cation process. Figure HJa) shows PL spectra at various 
pump powers for a device with A = 360 nm, now pro- 
duced in a low QD density region of the sample. Three 
isolated exciton lines arc observed on top of a broad back- 
ground near 942 nm. The sharp lines red-shift with in- 
creasing temperature (Fig. [3Jb)) with a dependence that 
can be fit to a model that predicts a red-shift of the InAs 
bandgap (Fig. |3tc)) 0, H, 0|- I n contrast, the broad 
background observed in Fig. EH a) shifts more slowly with 
temperature, and likely originates from out-coupling of 
broad QD multiexcitonic emission via the leaky cavity 
mode [ll| . This is reinforced by the observation, in 
Fig. EJb), that the sharp QD lines are maximized in the 
wavelength range 940 nm < A < 942 nm, when aligned 
to the broad cavity peak, and decrease when driven away 
from it. The slower cavity mode shift with temperature 
corresponds to a shift in refractive index [l2T |. 

Figure SJa) shows the detected PL as a function of av- 
erage pump power for the excitonic lines Xi and X2 and 
the cavity mode emission from Fig. [3ja). While Xi and 
X2 saturate at ss 20 nW, the cavity emission increases 
past this level. This further supports our assignments of 
QD transitions and cavity mode in the Fig. Ela) spec- 
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FIG. 4: (a) PL as a function of pump power for Xi, X2 and 
cavity emission from Fig. EI a), and two QDs in unpatterned 
GaAs. Error bars are 95% fit confidence intervals, (b) Solid: 
Xi lifetime trace with fit. Dotted: lifetime trace for QD em- 
bedded in a suspended GaAs waveguide. 

Improved photon extraction efficiency can potentially 
be achieved with a higher NA collection optic (50% in- 
crease for NA=0.7) and fabrication control [7j, while de- 
terministic QD spatial alignment [H], [H| can enhance 
both the efficiency and Purcell factor. Although single- 
photon emission from the bullseye is accompanied by un- 
desirable cavity emission, a few devices exhibited consid- 
erably less cavity mode feeding, albeit with lesser extrac- 
tion efficiencies. Since enhanced extraction efficiency is 
due to the directional far-field pattern, a trade-off may be 
achieved between Purcell enhancement and cavity feed- 
ing for reduced Q. It is also likely that quasi- resonant QD 
pumping will lead to reduced cavity feeding (l6j . These 
possibilities are under investigation. 

In summary, we developed a nanophotonic circular 
grating that provides ~ 10 % free space collection effi- 
ciencies for single InAs QD photons within a wavelength 
range of « 5 nm. Lifetime reduction of at least a fac- 
tor of four is achieved, which, taken together with the 
enhanced collection efficiency, indicates Purcell radiative 
rate enhancement. This structure allows for efficient and 
broadband spectroscopy of single QDs, and has potential 
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for use as a bright single-photon source. 
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SUPPORTING INFORMATION 



1. Design simulations 

The following simulation results illustrate the effects of varying trench depths on emission properties of the circular 
dielectric grating. Figure [if a) shows total emitted power P t ot as a function of wavelength for various trench depths, 
and Fig.[]Tb) shows the corresponding upwards (-P z +) and downwards (P z ~) extracted powers. 

Apparent in Fig. 0Ja) are a significant central wavelength shift and a strong radiative rate modification. Indeed, 
as shown in Fig. [TJc), the central wavelength blue shifts more than 40 nm for depths increasing from 95 nm to 190 
nm. The resonance full width at half maximum (FWHM), correspondingly, decreases, indicating an increase in field 
confinement and, consequently, the Purcell Factor (F p ), Fig. [He). The increased field confinement for deeper trenches 
is a consequence of better overlap of the guided field inside the slab and the etched region, which leads to increased 
guided wave reflectivity and reduces coupling to out-of-plane radiation. Figures QJb) and (f) show the effect of grating 
asymmetry on the ratio between upwards and downwards emitted powers. For a symmetric grating with d = t, 
emission in both directions is the same. Upwards emission is maximized for d/t « 0.8, for which P z+ /P z - = 2.19. 
This ratio increases to 5.8 for d/t = 0.7, however for a reduced F p . Clearly, a trade-off must be reached between 
Purcell enhancement and asymmetric emission. 



2. Effective mode volume calculation 

The bullseye cavity's effective mode volume T4ff=1.29(A c /n) 3 quoted in the text was obtained with the expression 



V eS = 



Je(r)|E(r)|Vr ^ 
max je(r) |E(r)| 2 |' 



(1) 



where E is the modal electric field and e(r) the medium permittivity. From V e s and the calculated cavity Q ~ 200, 
we can estimate the maximum Purcell enhacement factor F p [l[: 
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FIG. 1: (a) Total emitted power P tot as function of wavelength for various trench depths, (b) Vertically emitted power ±z 
direction (continuous:P z+ , upwards; dotted:P z _, downwards), as a function of wavelength, (c) Central wavelength, (d) full 
width at half maximum, (e) Purcell enhancement factor and (f) maximum upwards (P z +) and downwards (P z -) emitted power 
as functions of trench etch depth d. Phom is the homogeneous space electric dipole emitted power and t is the GaAs slab 
thickness. 



where n is the refractive index and A c is the cavity mode wavelength. The value determined through this calculation, 
F p = 11.8, corresponds well with the value F p = 11.0 that we determine from simulations of the enhanced radiated 
power for an electric dipole in the bullseye. The above values for F p assume perfect dipole orientation with respect 
to the cavity field and optimal dipole location within the field (i.e., in the bullseye center). 



3. Fabrication 



Based on simulation parameters, gratings were fabricated on a t=190 nm thick GaAs layer containing a single layer 
of self- assembled In As quantum dots on top of a 1 thick AlGaAs sacrificial layer. The epiwafer was grown with 
molecular beam epitaxy, and displayed a quantum dot density gradient from > 100 /im~ 2 to /im~ 2 along the (Oil) 
direction. Electron-beam lithography was used to define the patterns, and a single, timed, inductively-coupled plasma 
reactive ion etch (ICP-RIE) step with an Ar/Cl2 chemistry transferred the gratings into the GaAs. This step was 
optimized so that the GaAs would be partially etched to a desired depth in the grating region (see Fig. 1(b)), and 
fully etched over the large, curved rectangles just outside the grating region seen in Fig. 1(a). These open areas were 
included to give access to the AlGaAs sacrificial layer, which was etched with Hydrofluoric acid in a final step. Devices 
with varying quantum dot (QD) densities were produced by fabricating the devices along the QD density gradient of 
the wafer. 



4. Photoluminescence spectra and collection efficiency 



The spectra shown in Fig. 3 were obtained with a grating spectrometer and a Si charge-coupled device (CCD). To 
obtain the PL intensity from the excitonic lines Xi and X2 in Fig. 3(a) without contributions from the broad cavity 
background, Lorentzians were fitted to the corresponding peaks. Emitted photon rates plotted in Fig. 4(a) for Xi 
and X2 and for unpatterned GaAs QD lines correspond to the Lorentzian areas. The integrated cavity photon rates 
in Fig. 4(b) were obtained by integrating the spectra between 930 nm and 955 nm and subtracting rates from the 
two excitonic lines. 

To estimate the collection efficiency, we convert detected CCD counts to photon counts into our NA=0.42 objective. 
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The two are related by a conversion factor that is equal to the product of our detection efficiency £, the transmission 
through the PL setup T pat h, and the transmission through the cryostat windows T w ; n( j ows . £ includes the in-coupling 
efficiency into the spectrometer, the spectrometer's grating efficiency, and the CCD's quantum efficiency, and is 
determined by sending a reference laser with known power and wavelength into the spectrometer. In particular, we 
attenuate a 102 nW laser at 960 nm by 50 dB and acquire a spectrum with a I s integration time. Integrating the laser 
spectrum yielded a count rate of 6.46 x I0 4 s _1 , which, when compared to the photon rate just before the spectrometer 
(4.93 x 10 6 s _1 ), gives a factor of w 77 photons per count. Because the QD emission wavelength (Ax x = 941 nm) 
differs from the calibration wavelength of 960 nm, we multiply this conversion factor by 0.78, which is the ratio of 
the manufacturer-specified CCD quantum efficiencies at these two wavelengths. We therefore get an overall detection 
efficiency £ = 0.0167 (= 60 photons per count). 

The transmission through the optical path (including the collection objective) was T patn = 0.156, and was measured 
by launching a fiber-coupled laser of known power and wavelength through the optical setup. Note that, since the 
fiber NA=0.13 is less than the objective's NA, all of the power emitted from the fiber is collected, so the measured 
transmission includes only transmission losses through the objective and routing optics. T w i n dows ~ 0.87, and includes 
transmission through the radiation shield and outer cryostat windows. 

The QD line Xi in Fig. 4(a) yields a saturated CCD count rate of R Xl = 1.14 x 10 4 s _1 ±0.16 x 10 4 s~\ where the 
uncertainty is a 95 % fit confidence interval, due to spectrometer resolution and detection noise. Assuming 100 % QD 
quantum efficiency, the rate of single photons emitted by the saturated QD is equal to the pulsed pump excitation 
rate, like R cx = 50 MHz. The collection efficiency was then calculated as 

_ ^Xj ( o\ 

s ' J path ' windows 

Substituting all these values yields a collection efficiency 77 = 10.1 % ± 1.4 %. 



5. Quantum dot temperature dependence 

The temperature dependence of the sharp features X 1; X 2 and X 3 in Figs. 3(a) and (b) was fitted with the 
Bose-Einstein expression 



E res {T) = E res (T = 0) - Shu 



coth ——= 
\2k B Tj 



(4) 



This expression models the evolution of the semiconductor bandgap energy with temperature due to electron-phonon 
interaction, assuming no phonon dispersion, and has been successfully applied towards excitonic transitions in epitax- 
ially grown quantum dots [HQ. I n Eq. ((4]), T is the sample temperature, E res (T) is the excitonic resonance energy, 
huj is the phonon energy, and S is a dimensionless coupling constant. The fits shown in Fig. 3(c) were obtained with 
the parameters in Table HI 

The phonon energies Hoj in Table U are between the bulk GaAs transverse acoustic phonon energies at the X and 
L points, HojtaOQ = 7.7 mcV and HutaO^) — 9.8 meV, and the coupling coefficients S are compatible with those 
reported in refs. and [3- This indicates that the sharp spectral lines correspond to excitonic QD transitions. 





Xi 


x 2 


x 3 


E res {T = 0) (eV) 


1.3186 ±0.7 x 1(T 4 


1.3173 ± 0.4 x 10~ 4 


1.3161 ±0.3 x 10" 4 


hu) (meV) 


9.1601 ±0.8632 


9.1373 ±0.4890 


8.7596 ± 0.4536 


S 


0.9093 ±0.1005 


0.9276 ± 0.0581 


0.8662 ± 0.0581 



TABLE I: Fitting parameters for temperature dependence of QD lines Xi, X2 and X3 with eq.©. Errors are 95 % fit confidence 
intervals (two standard deviations). 
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